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Abstract 

We present a preliminary measurement of the branching fraction for B — > peu and of the 
CKM matrix element \V U b\ using approximately 55 million BB meson pairs collected with the 
BABAR detector. Using isospin relations for several modes we find 

B(B° -» p~ e+ v) = (3.39 ± 0.44 ± 0.52 ± 0.60) x 10~ 4 
\V ub \ = (3.69 ± 0.23 ± 0.27 x 10" 3 . 

The quoted errors are statistical, systematic, and theoretical respectively. These results are 
obtained by using five different form-factor calculations. 
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1 Introduction 



Exclusive b — > u£f decays can be used to determine the modulus of V u b, one of the smallest 
and least well known CKM matrix elements. Compared to the determination with inclusive 
decays, the extra kinematical constraints allow access to a larger part of the lepton-momentum 
spectrum, resulting in smaller extrapolation uncertainties. Experimentally, the main difficulty 
for the observation of b — * utv signal events is the large background from b —* civ events. 
Because IV^/V^I ~ 0.1, the branching fractions of the exclusive b —* ulv decays (~ 10 ) are 
small compared to those of the charmed semileptonic decays, which are of the order of some 
percent. 

To extract the branching fraction B — > pev and \V u h\ requires use of hadronic form- factors 
which have to be obtained from theory. In this analysis we use five different form factor cal- 
culations: the two quark models ISGW2 [GJ and Beyer /Melikhov the lattice calculation by 
the UKQCD group ||, a model based on light cone sum rules (LCSR §), and a model based 
on heavy quark and SU(3) symmetries (Ligeti/Wise ||). The UKQCD and LCSR calculations 
directly use QCD, whereas the quark models are more phenomenological. All calculations use a 
particular value of q 2 = (p e + p u ) 2 as a normalization point. Usually, is used as this is the 
point where the hadronic system is least disturbed. The LCSR result however is normalized at 
a lower value of q 2 . 

2 The BABAR Detector and Data Set 

The data used in this analysis were collected with the BABAR detector || at the PEP-II e + e~ 
storage ring The integrated luminosity of the sample is 50.5 fb _1 taken at the T(AS) mass 
( "on-resonance" ) , corresponding to 55.2 million BB meson pairs. An additional 8.7 fb -1 of data 
were taken 40MeV below the T(4S) resonance ("off-resonance"). 

PEP-II is an e + e~ collider operated with asymmetric beam energies, producing a boosted 
(/?7 = 0.55) T(4S) along the collision axis. BABAR is a solenoidal detector optimized for the 
asymmetric beam configuration at PEP-II. Charged particle (track) momenta are measured in a 
tracking system consisting of a 5-layer, double-sided, silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) filled with a gas mixture of helium and isobutane, both operating within 
a 1.5 T superconducting solenoidal magnet. Photon candidates are selected as local maxima of 
deposited energy in an electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals 
arranged in barrel and forward endcap subdetectors. Particle identification is performed by 
combining information from ionization measurements (dE/dx) in the SVT and DCH, and the 
Cherenkov angle 9 C measured by a detector of internally reflected Cherenkov light (DIRC). The 
DIRC system is a unique type of Cherenkov detector that relies on total internal reflection within 
the radiating volumes (quartz bars) to deliver the Cherenkov light outside the tracking and 
magnetic volumes, where the Cherenkov ring is imaged by an array of ~ 11000 photomultiplier 
tubes. The detector is surrounded by an instrumented flux-return (IFR). 

3 Event Selection 

In this section, we describe the selection of the exclusive semileptonic decays B + — > p°e + v, 
B° -» p~e + v, B + -> we + i/, B+ -» 7r°e + v, and B° -» Tr'e+v (with p° -> vr+vr-, p ± -> 7r°vr ± 
and oj — ► 7r°7r + 7r~). The charge conjugate decays are implied throughout. Our analysis strategy 
is similar to one used by CLEO ||. The analysis is optimized for B — > pev decays; the it and 
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ijj modes are included because of the crossfeeds into the p modes. Isospin and quark model 
relations are used to effectively measure only two branching fractions, one for B — > peu and 
one for B — ► irev. This is described in section [|. This analysis uses only electrons and not 
muons because the background contribution from fake leptons is much lower in the case of 
electrons. We reconstruct three kinematic variables which are used in the fit to extract the 
signal yields. These are the electron energy E^^, the invariant hadronic mass M 7rn ^ (for the 
p and co modes), and the difference between the reconstructed and expected B meson energy 
(AE = -Ehadron + E t + |p m iss|c - E beam ) in the center-of-mass (CM) system. 

Two electron energy regions are considered: 2.0 < Eg* < 2.3 GeV (LOLEP), and 2.3 < 
Eg* < 2.7 GeV (HILEP). The HILEP region is most sensitive to the signal because the b — > ceu 
events are almost completely suppressed; the largest background source here is from continuum 
e + e~ — ► qq events. Real data taken below the T(45) mass, which includes e + e~ — ► e + e~( r y) 
events, is used for the continuum subtraction. In the LOLEP region, b — > c decays dominate 
and provide the normalization of the background at higher electron energies. 

Hadronic events are selected based on track multiplicity and event topology. The tracks 
must have at least 12 hits in the drift chamber. We also require that the impact parameter 
of the track along and transverse to the beam direction must be less than 3 cm and 1 cm, 
respectively. In addition, the transverse momentum must be greater than 0.1 GeV/c. Clusters 
in the electromagnetic calorimeter of BABAR that are not associated to any tracks must have 
an energy greater than 30 MeV to be considered as photons. In addition, the lateral moment 
of the shower energy distribution |l0| must be smaller than 0.8. We select events with either 
-^tracks > 5 or (-/V trac k s > 4 and iV p hotons > 5). The B mesons are produced nearly at rest so 
their decay products are distributed roughly uniformly in solid angle. In contrast, continuum 
events have a much more collimated (jet-like) event topology. They are suppressed by requiring 
the ratio of Fox- Wolfram moments H2/ Hq || to be less than 0.4. This requirement suppresses 
55% of the e + e~ — ► qq background and non-hadronic events, with a signal efficiency of 85%. In 
addition a neural net is used for further suppression of continuum events, as described below. 

To identify electrons, we use a likelihood estimator, which uses information from several 
BABAR sub-detectors. The primary information is the ratio of the calorimeter energy to the 
track momentum. We require that the direction of the electron momentum is within the good 
calorimeter acceptance —0.72 < cos6> e) i a b < 0.92. The efficiency of this selector is around 90%, 
with a pion misidentification rate of less than 0.1%. We also reject electrons from J/i/j decays 
(requiring two electrons identified with the likelihood electron selector and 3.00 < M e + e - < 
3.14 GeV/c 2 ) and from photon conversions (M e + e - < 30 MeV/c 2 ). 

To reconstruct the neutral p meson we combine two oppositely charged tracks, and for the 
case of the charged p a track and a 7r°. The 7r° mesons are reconstructed from two photons with 
an invariant mass 120 < M 77 < 145 MeV/c 2 corresponding to ±2 a from the nominal ir° mass 
on average. To suppress combinatorial backgrounds we require that the pion with the highest 
momentum must have p^ M > 400 MeV/c and the other pion must satisfy p^ M > 200 MeV/c. 
For the u, we combine two oppositely charged tracks with a tt°. The ui invariant mass is measured 
within ±80 MeV/c 2 of the nominal lo mass [11]. This includes a side band region below and above 
the nominal lo mass. To suppress combinatorial backgrounds we require > 100 MeV/c for 
each of the three pions. The charged tracks used to reconstruct the rho, omega, or pi^ mesons 
must not have been identified as kaons. 

In the following discussion all variables are taken in the center-of-mass frame. The neutrino 
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is reconstructed from the missing momentum: 

Pmiss = - ^2 Pi- ^2 Pi , (!) 

tracks photons 

where the sums run over all reconstructed and accepted tracks and photons in the event. We 
then take (E u ,p u ) = (\p m i SB \c, 

Pmiss ) ■ 

A 5-meson decay consistent with the signal modes is reconstructed using the constraints 
Eb = -Ebeam and (pb — Py) 2 = 0, where Y is the (p, u>, tt) + e system. A useful quantity for 
testing consistency is the angle between the B momentum direction and that of the reconstructed 
Y system, see Fig. |], 

_ 2E B E Y -{Ml + M Y )c± 

COS 9 BY - ||^ I o • V') 

2\pb\\py\c z 

Background tends to have non-physical values of cos##y. For the signal, small extensions of 
| cos#ey| > 1 are allowed because of detector resolution. We therefore require 

\cos6 B y\ < 1-1 • (3) 

The efficiency of this requirement is almost 100% for the signal and it rejects more than 60% of 
the b — > ceu and 80% of the continuum background. 

We also compare the direction of the missing momentum p^iss with that of the neutrino 
momentum inferred from p u = Pb — Py- The latter is known to within an azimuthal ambiguity 
about the ps direction because the magnitude, but not the direction, of the B meson momentum 
is known. We use the smallest possible angle A6> m i n between the two directions, which is obtained 
when the momenta of py, p u , arid p m i ss are in the same plane, see Fig. [l|. We use the requirement 

0.8 < cosA0 min < 1.0 . (4) 

This has been optimized using a Monte Carlo simulation, in such a way as to minimize the 
relative error on the measured branching fraction. 
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In addition, we require | cos# m i ss [ < 0.9 where 6 m i ss is the angle between p m i ss and the beam 
axis. This cut rejects events with missing high momentum particles close to the beam axis. 

The continuum e + e~ — ► qq, where q = u,d, s, c, is an important background at high electron 
energies, where we are most sensitive to the signal. To reject these events, we use a neural net 
with 14 event shape variables such as the track and cluster energies in nine cones around the 
electron-momentum axis. The optimized cut on the neural net output suppresses more than 
90% of the continuum, after all other requirements have been applied, in the HILEP region. 
The selection efficiency on the signal is 60%. 

After all the above criteria, there can still be several candidates per event. This follows from 
the large width of the p and also because we are reconstructing five different decay modes. To 
avoid statistical difficulties related to large numbers of combinations, we choose one candidate per 
event, namely the one with a reconstructed total momentum closest to the B meson momentum: 

IPhadron + Pe + Pmiss | closest to \p B \- (5) 

The efficiency of this selection for the signal is close to 85% . 

At each step of the selection procedure the data distributions agree well with their Monte 
Carlo simulation. 

The total efficiencies in the HILEP region are 4.21% for the B+ -» p°e+is mode and 3.31% 
for the B° — > p~~e+v mode. These efficiencies are determined using the ISGW2 form-factors. 
They are defined here as the number of HILEP signal events that pass all selection criteria and 
are reconstructed in the specified channel divided by the number of all generated events (of all 
electron momenta) for this channel. The signal and crossfeed efficiencies for all channels are 
listed in Table [|. 

Table 1: Selection efficiencies for the modes B+ — > p°e+u, B° — > p~e+u, B+ — > oje+v, B+ — > 
and B — > ir e+v in the HILEP region. These efficiencies have been determined using 
simulated data (ISGW2 model). They are defined as the number of selected events passing all 
cuts divided by the number of all generated events (in the full lepton energy range) for the 
specified channel. 



Reconstruction Efficiency (%) 



Generated 


B+ p^e+u 


B [) -> p~e + v 


B + -> Lue+v 


B+ -» 7r u e+z/ 


B [) -> TT'e+v 


B+ -> p^e+v 


4.21 


0.99 


0.30 


0.58 


1.35 


B° -> p-e+v 


1.04 


3.31 


0.31 


1.08 


1.32 


B+ -» ue+v 


1.97 


1.55 


1.57 


0.94 


0.88 


B+ _» 7r e+zy 


0.23 


0.41 


0.07 


1.28 


0.40 


B° -> Ti-e+u 


0.42 


0.47 


0.05 


0.33 


1.63 



4 Fit Method 

We have performed a binned maximum-likelihood fit of the two-dimensional distribution (M^M > 
AE) simultaneously in the two electron energy ranges (LOLEP, HILEP) and the decay modes 
B+ — > p°e+u, B° — > p~e+v, and B+ — > uoe+v. For the B — > pev modes, the data are divided 
into 10 x 10 bins over the {M n7r , AE) region 0.25 < < 2.00 GeV/c 2 and \AE\ < 2 GeV. 
The bin size for the fit is thus 175 MeV/c 2 in and 400 MeV in AE. For the oj channel, we 
use 5 bins in the range 782±80 MeV/c 2 and 10 bins in \AE\ < 2 GeV. The modes B+ -> iPe+v 
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and B° —* ir e + v are also included to model the crossfeeds into the other signal channels; for 
these modes only AE is used clS £L fit variable. 

Our fit includes contributions from the signal modes, other b — » uef decays, 6 — » cev decays, 
continuum, and a contribution from misidentified electrons. For the signal and backgrounds 
coming from other b — > uef and b — > cef decays, Monte Carlo simulation provides the shapes of 
the distributions. The decays B — » have been simulated using heavy quark effective theory 
(HQET Ifl^l ). The modes B — » D*ir£v are simulated according to the Goity- Roberts model [jl3|| . 
Resonant 6 — > u£i^ downfeed modes are implemented according to the ISGW2 model. Non- 



T{B° - 


■* p e + v) 


= 2T(B + 


-> p°e+v 


T{B + - 


-+ p°e + v) 


= T{B + - 


-> toe + v), 


r(B° - 


* ir~e + v) 


= 2T(B + 





resonant b — > u&/ modes are implemented according to a model by Fazio and Neubert [14]. 

Isospin and quark model relations are used to constrain the relative normalizations of B° — > 
p~e + v, B + —> p°e + v, and -B + — > we + y and therefore to reduce the number of free fit parameters: 

(6) 
(7) 
(8) 

Isospin breaking effects are discussed in fl5|| and [|lfj]. We assume that the isospin relations in 
Eqs. |TJ and |^ are broken by not more than 3%. This would have a negligible effect on our result 
and therefore we do not include a corresponding systematic error. The isospin relations were 
tested experimentally. We find that the isospin relations Eqs. |fj and fj] are consistent with the 
data within 1.3 a and 1.7 a. 

We use the following 9 free parameters for the fit: 

• B(B° — ► p~e + v) (1 parameter); 

• B(B° — > Tt~e + v) (1 parameter); 

• the scale factors of the b — > uev background in each electron energy bin (2 parameters), 
that give the overall normalization of all b — > uev modes that are not signal modes, relative 
to that expected from the Monte Carlo simulation; 

• the scale factors, one for each mode, that give the overall normalization of the b — > cev 
background relative to that expected from the Monte Carlo simulation (5 parameters). 

The maximum-likelihood fit method used in this analysis has been described in [17]. The 
fit takes into account the statistical fluctuations not only of the on-resonance and off-resonance 
data but also those of the Monte Carlo contributions. We have performed a toy Monte Carlo 
check to verify the stability of the fit method and to check the statistical error returned by the 
fit. 

Whereas the signal modes B + — > p°e + v, B° — > p~e + v, and B + — > uje + v are simulated with 
five different form- factors, all other b — > uev modes (downfeed background) are simulated using 
the ISGW2 form- factor only. 



5 Fit Results 

The signal yields extracted from the binned maximum-likelihood fit in the HILEP region are 
324 ± 40 B+ -> p°e + v events and 510 ± 63 B° -> p~e + v events, based on the ISGW2 
calculation. The composition of events for the B° — > p~ev and B + — > p°e + v channel is shown in 
Table ||. The isospin-constrained results for the five different form-factors are shown in Fig. [2|. A 
X 2 test has been performed to check the quality of the fit. Bins in sparsely populated regions have 
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Table 2: Summary of data yields for the B° — > p e + u and B + — > p G e + u modes with electron 
energies between 2.3 and 2.7 GeV (HILEP), and between 2.0 and 2.3 GeV (LOLEP). The yields 
presented in this table were obtained using the ISGW2 form- factor. The downfeed background 
includes all B — * X u ev modes except for p, u, and ir. The crossfeed signal contribution corre- 
sponds to events from the other signal modes with p , u, or n and is constrained to the signal 
in the fit. All errors are statistical only. 







-» p e + v 


B+ - 


-> p°e + u 




HILEP 


LOLEP 


HILEP 


LOLEP 


On-resonance yield 


2302 


39349 


2213 


40155 


Direct signal 


510 ± 63 


718 ± 89 


324 ± 40 


440 ± 55 


Crossfeed signal 


262 ± 32 


538 ± 73 


363 ± 42 


725 ± 86 


Downfeed 


203 ± 55 


2278 ± 403 


226 ± 92 


2435 ± 430 


b — ► ceu 


414 ± 5 


33859 ± 438 


367 ± 5 


34366 ± 458 


e + e~ — > qq 


917 ± 73 


1928 ± 106 


912 ± 73 


2063 ± 110 


Fake electrons 


12 ± 3 


80 ± 9 


18 ± 4 


76 ± 9 



been combined before the x 2 calculation. For ISGW2, we obtain x 2 = 91 f° r -Wdof = 93, which 
corresponds to a p-value of 0.52, and similarly good fit quality for the other four form-factor 
calculations. 

The five fit parameters describing the b — > c backgrounds agree with the Monte Carlo expec- 
tations within 9% on average. The two parameters describing the b — > u downfeed background 
in LOLEP and HILEP agree to better than 1.5 a and 1.2 a. The fit result for the n modes is 
B{B° -» 7r-e + u) = (1.87 ± 0.56) x 10~ 4 for the ISGW2 calculation. 

The projections of the ISGW2 fit result for the two electron energy bins after subtraction 
of the continuum contribution are shown in Figs. ^ and ||. Good agreement between the data 
and fit result is seen in each of these figures. The fits for the other four form-factor calculations 
show similar agreement. 



6 Systematic Errors 

The summary of all systematic errors on the branching fraction that have been considered is 
shown in Table ||. The total systematic error is taken as the quadratic sum of all individual 
errors. The fluctuations due to finite Monte Carlo statistics are included in the statistical error 
and not in the systematic error. 

The largest single systematic error comes from the uncertainty in the shape of the downfeed 
background. We use the ISGW2 model to describe resonant downfeed modes, and a model 



by Fazio and Neubert [14| for non-resonant modes. The fraction of non-resonant events in 
the downfeed background is varied from 0% to 68% to estimate this systematic uncertainty. 
The composition of the resonant b — > u downfeed component has been varied by changing the 
branching fraction for individual resonances by ±50%, and keeping the total rate constant. 

We have also varied the most important selection requirements of this analysis within a 
reasonable range and have changed our fit method (fitting with only four channels, or without the 
LOLEP region, or with different binnings). Most variations seen are within la of the expected 
statistical variation, some are close to 2a. To be conservative we assign a systematic error 
corresponding to half the largest variations seen. This corresponds to the last two systematic 
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ISGW2: 

2.79 ± 0.35 ± 0.43 

UKQCD: 

3.45 ± 0.44 ± 0.53 

LCSR: 

3.99 ± 0.52 ± 0.62 

Beyer/Melikhov: 
3.77 ± 0.48 ± 0.58 

Ligeti/Wise: 
2.93 ± 0.41 ± 0.45 



Combined: 

3.39 ± 0.44 ± 0.52 ± 0.60 

2 2.5 3 3.5 4 4.5 5 
Br(B°H> p"e + v)x 10 4 

Figure 2: The B° — > p~e + v branching fraction results using the ISGW2, UKQCD, LCSR, 
Beyer/Melikhov, and Ligeti/Wise form- factors. The errors shown are statistical, systematic, 
and (only in case of the combined result) theoretical, successively combined in quadrature. The 
combined central value is determined by taking the unweighted mean of all form-factor results. 
The statistical and systematic errors of the combined result are determined by taking the means 
of the relative errors of each individual result, and its theoretical error is taken to be one half 
of the full spread of the results. 

errors quoted in Table ||. 

7 Extraction of \V u b\ 

The CKM matrix element \V u b\ can be obtained from the branching fraction B(B° — > p~e + v) 
using 

K^J B(B "r^\ ( 9) 

V 1 thyT B° 

where Tth y is the predicted form-factor normalization. Values of r\hy and theoretical errors for 
each form-factor calculation are given in Table The calculations quote errors on Tthy between 
15% and 50%. We use t#o = 1.548 ± 0.032 ps [O], and the branching fractions are taken 
separately for each form-factor as listed in Fig. 0. The combined central value is determined by 
taking the weighted mean of all form- factor results. The statistical and systematic errors of the 
final combined result are determined by taking the mean of the relative errors of each individual 




13 



9000 




0.4 0.6 0.8 1 1.2 1.4 1.6 1.8,2 
M had [GeV/c ] 



1.2 1.4 1.6 1.8 ,2 

M had t GeV/c ] 



(a) Mnn (LOLEP) 



(b) (HILEP) 




0.5 1 1.5 2 
A E [GeV] 



0.5 1 1.5 2 
A E [GeV] 



(c) AE (LOLEP) 



(d) AE (HILEP) 



Figure 3: Continuum-subtracted projections of the ISGW2 fit result for the B + — > p°e + u 
channels in the LOLEP and HILEP electron energy regions; the contributions are the direct 
and crossfeed components of the signal (unhatched region, above and below the dashed line, 
respectively); the background from b — > uev other than B — > pev and B — > ujev modes (double- 
hatched region); the background from b — ► ce^ and other backgrounds (single-hatched region). 
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Figure 4: Continuum-subtracted projections of the ISGW2 fit result for the B° — > p~e + v 
channels in the LOLEP and HILEP electron energy regions; the contributions are the direct 
and crossfeed components of the signal (unhatched region, above and below the dashed line, 
respectively); the background from b — > uev other than B — > pev and B — > wef modes (double- 
hatched region); the background from b — ► cef and other backgrounds (single-hatched region). 
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Table 3: Summary of all contributions to the systematic error on the branching fraction for 
B — > pev. 



Error contribution 




Irackmg x^mciency 


±5 


Tracking Resolution 


±1 


Photon/ 7T U bmciency 


±5 


Photon/ 7r Energy bcale 


±3 


b — > c Background Composition 


+1.4 
-1.7 


Resonant b — > it Background Composition 


+6 
-4 


i> on-xvesoiiaiiL c* — > u j3dCKgiou.no 


-l-Q 

zty 


B Lifetime 


±1 


B Counting 


±1.6 


Fake Electrons 


< ±1 


Electron Id 


±2 


/±//oo 


< ±1 


Data Selection 


±6 


Fit Method 


+4 

-6 


Total Systematic Error 


±15.5 



Table 4: r t h y predicted by various form-factor calculations. 



Form-factor f t h y (ps x ) Estimated error on fth y (%) 

ISGW2 142 

LCSR 16.9 

UKQCD 16.5 

Beyer /Melikhov 16.0 

Ligeti/Wise 19.4 



result. The theoretical error is taken to be one half of the full spread of all fit results (including 
theoretical errors). The results for each form- factor and the combined result is shown in Fig. |5|. 

A comparison of our preliminary result with inclusive and exclusive measurements from 
CLEO and the inclusive measurement from LEP is shown in Fig. ^. Two exclusive results from 
CLEO are quoted. The first result is obtained from an analysis very similar to the analysis 
presented here || , the second result is an average of their first result and a separate analysis || 
[l8| . Our result is compatible with all other measurements within errors and lies between the 
CLEO and LEP results. 



±50 
±32 

+21 
-14 

±15 
±29 
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3.5 



4 4.5 
IV J x 10- 



ISGW2: 

3.56 ± 0.22 ± 0.27 

UKQCD: 

3.68 ± 0.23 ± 0.27 

LCSR: 

3.91 ±0.25 ±0.29 

Beyer/Melikhov: 
3.90 ± 0.24 ± 0.29 

Ligeti/Wise: 
3.12 ±0.21 ±0.23 



Combined: 



3.69 ± 0.23 ± 0.27 



+0.80 
-1.04 

+0.37 
-0.27 

+0.58 
-0.68 

+0.28 
-0.30 

+0.42 
-0.49 



+0.40 
-0.59 



Figure 5: \V ub \ determined using the ISGW2, UKQCD, LCSR, Beyer/Melikhov, and 
Ligeti/Wise form- factors. The results for each form- factor are drawn with theoretical error 
bars only. The combined central value is determined by taking the mean of all form-factor re- 
sults, weighted by their individual theoretical error. The theoretical error of the combined result 
is taken to be one half of the full spread of results, including the errors. The combined result is 
also drawn with statistical, systematic, and theoretical errors successively added in quadrature. 
The statistical and systematic errors of the combined result are determined by taking the mean 
of the relative errors of each individual result. In addition we give, on the right side of the figure, 
all five results with their statistical, systematic, and theoretical errors. 
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2.5 3 3.5 4 4.5 



BABAR exclusive, this result 
3.69 ± 0.23 ± 0.27 +° 5 4 9 ° 



CLEO exclusive I 
3.23±0.24 +° 2 2 6 3 ± 0.58 



CLEO exclusive II 
3.25 10.14 ^±0.55 



CLEO inclusive: 

4.08 ± 0.63 (stat+syst+theo) 



LEP inclusive: 

4.09 ± 0.68 (stat+syst+theo) 



IV J x 10 : 



Figure 6: Comparison with results from other experiments. The CLEO exclusive I result @ is 
obtained from an analysis very similar to the analysis presented here, their exclusive II result 
fD| is an average of the exclusive I result and a separate CLEO result. 
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